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Abstract. We consider two hidimensional classical Ising mod- 
els, coupled by a weak interaction bilinear in the energy densities 
of the two systems; the model contains, as limiting cases, the 
Ashkin-Teller and the Eight-vertex models for certain values of 
their parameters. We write the energy-energy correlations and the 
specific heat as Grassman integrals formally describing Dirac l-hl 
dimensional interacting massive fermions on a lattice, and an 
expansion based on Renormalization Group is written for them, 
convergent up to temperatures very close to the critical tempera- 
ture for small coupling. The asymptotic behaviour is determined 
by critical indices which are continuous functions of the coupling. 
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1. Introduction and main results 



1.1 Motivations. It is well known that the partition function and many correla- 
tions of the d = 2 classical Ising model can be written as Grassman integrals for- 
mally describing a system of d = free fermions, see [H] , [SML] , [K] , [MW] , [S] , [ID] . 
This mapping of the Ising model in terms of a fermionic theory is quite useful; as 
the corresponding action is quadratic in the Grassman variables, it can be diago- 
nalized by a Bogoliubov transformation so that the partition function and many 
observables can be exactly computed. However, even a slight perturbation of the 
Ising model makes the corresponding fermionic action non-quadratic. The rela- 
tionship between spin Ising-like models and interacting fermionic models is deeply 
investigated in the literature; it is for instance claimed that models of coupled 
Ising models with an interaction quartic in the spins, like the Eight-vertex model, 
are in the same class of universality of models of interacting d = 1 + 1 fermions in 
the continuum, like the massive Thirring model or the massive Luttinger model, 
see for instance [N], [LP], [PB]. This relationship is however found under many 
approximations. In this paper we consider two Ising models coupled by a general 
interaction bilinear in the spin densities; the model contains as particular cases 
the Ashkin- Teller and Eight-vertex model for certain values of their parameters. 
We clarify the relation between this general spin model and systems of interacting 
fermions by writing the partition function and some correlations as Grassman in- 
tegrals with a formal action resembling but not coinciding with the action of the 
massive Thirring or Luttinger model; the main difference (in addition to the fact 
that the formal action describes fermions in a lattice, not in the continuum) is 
that it does not verify some special symmetries (like invariance under local Gauge 
transformations) which are generally verified in models of interacting fermions 
and which play an important role in their analysis. The main interest in this 
representation in terms of Grassman integrals is that we can apply the Renormal- 
ization Group techniques developed for interacting fermions (see for instance [BG] 
or [GM] ) for writing a convergent perturbative expansion for the partition function 
and some correlations of the spin model we are considering. A naive power series 
expansion in the coupling is not expected to be convergent close to the critical 
temperature (the free and interacting model are not analytically close) and one 
needs a resummed or renormalized expansion which is provided by Renormalization 
Group. 

1.2 Spin systems with quartic interactions. We consider two Ising models coupled 
via a four spin interaction bilinear in the energy densities of the two sublattices. 
Given Am G a square lattice with side M and periodic boundary condition, we 
call X = (x, xq) a site of Am- If (T'x^'* = ±1 and = ±1, we write the following 
Hamiltonian 

if^(a(i), a(2)) = Hi{a^^^) + Hi{a^^^) + V{a^^\a^^^) (1.1) 
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where, if a = 1,2 

M 

Hiia(-^) = - E [ji'^VilaS,., + (1.2) 

a;,a;o=l 

M 

F((J*^^\ (7*^^^) = — ^ {•^^['^i^io^a;+l,xo'^a;^Ko'^a;+l,xo + '^a;!a;o'^a;,a;o + l^a;!a;o^a;,Ko+l] 

+AfofrT(l) a^^^ rr^^) (2) (i) (l) (2) (2) . . o^ 

^'^"[^x,xo"x+l,xo^x,xo^x,xo+l ^ ^x,xo^x,xo+'i-^x-l,xo+l^x,xo + li y^-"^) 

L'^a;,a;o'^a;+l,a;o'^a;,xo'^a;,xo + l ^ a;,a;o"a;+l,xo"a!+l,a;o-l'^a;+l,a;oJ 

a=l,2 

The above Hamiltonian can be a model for a number of physical problems; for 
instance two magnetic spin planes coupled by an interplanar interaction quartic 
in the spins (in this interpretation a^^^ and cr^^^ are the spin of the first or second 
plane). Moreover it contains, as limiting cases, two well known models, the Ashkin- 
Teller and the Eight-vertex model. 

1.3 The Ashkin-Teller model. Ashkin and Teller introduced their model as a 
generalization of the Ising model to a four state system. Each site x is occupied 
by one of the four kinds of atom, said A,B,C,D, and two neighboring atoms interact 
with an energy Sq for AA, BB, CC\ DD; ei for AB, CD, 62 for AC, BD, for 
AD,BC. It is quite easy to express this model in terms of Ising spins, see [F]. 
One associates to each site of the lattice two spins and crx \ see Fig.l; then 
o'x\o-x^ = (+,+) if there is an A atom at site x, (+, — ) if there is an atom 
B, (— ,+) if there is an atom C and (— ,— ) if there is an atom D. Then the 
hamiltonian is given by 

iy^(aW,^(2)) = iy,(^(i)) + if,(^(2)) 

M 

-J" V [(7^2) (2) ^(1)^(1) +^(2)^(2) (jW 1_J^ (14) 

X,Xo = l 

where Hi{a^'^^) is given by (1.2) with J^"^ = J^"^ = and 

- J^i) = (£0 + £i - £2 - e3)/4 - J(') = (£o + £2-£3-ei)/4 (1.5) 

-J" = (£o + £3 - £i - £2)/4 - Jo = {eo + £i + £2 + £3)/4 

The model is not soluble, except in some particular case, and we refer to [Ba] 
chapt.l2 and references therein for a review of the main results about it; it is of 
the form (1.1) with b — c — 0. 

1.4 The Eight-vertex model. In such model one associates to each site of a square 
lattice a vertex with four arrows with different orientations, see [Ba] chapt. 10. 
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Among the sixteen possibilities, only eight kind of vertices are allowed, and an 
energy Sj, z = 1, .., 8, is associated to each of them. It was solved in [Ba] if £i = £2, 
£3 = £4, £5 = £6, £7 = £8 and for suitable values of the Si it reduces to the ice- 
model, solved in [Li]. Such models are called vertex models and were introduced 
to describe crystal with hydrogen bonding. As explained in [Ba] or [W] they can 
be written in terms of the following Hamiltonian, see Fig. 2 

HA{cr^^\ O"^^^) = - (Tx,xo+l(^x+l,xo + J'(^x,xo(^x+l,XQ + l 

x,xo 

+ J" 

^X,Xo^X+l,Xo + 'i-^X,Xo + 'i-'-^X+l,Xo ] (1-6) 

with 

£1 = £2 = - J - J' - J" £3 = £4 = J + J' - J" (1.7) 

£5 = £6 = J' - J + J" £7 = £8 = - J' + J + J" 

If J" = the above hamiltonian decouples in the hamiltonian of two Ising models, 
involving spins located in two sublattices, one with x + xq equal to an even number 
and the other to an odd number. Each site of the sublattices is in the center of 
the unit square of the other one. Calling a^^^ and cx*^^-' the spin located in the first 
or the second sublattice, relabeling the spins and performing a rotation of ^ the 
hamiltonian can be written as, see Fig. 3 

iy^(^U),^(2)) = iy,(^(i))+iy,(^(2)) 

M 

- I" V [^7^2) (2) (1) (1) (1) (1) (2) (2) , . 

X,Xo='^ 

where 

M 

X,Xo = l 

which is again of the form (1.1) with a = c = (of course this identification is 
exact if boundary conditions in the Eight- vertex model are chosen properly). 

1.5 Main results. We will consider in this paper the case 

jf ^ = 4^ = Jf ) = ) = J (1.9) 

If A = the problem reduces to the Ising model whose critical temperature is 
given by tanh /9c Jo = — 1 = tc where J = /9Jo and f3~^ is the temperature. The 
interaction changes in general the value of the critical temperature by terms 0(A), 
and it is instead technically convenient to fix its value still in correspondence of 
tc = \f2 — 1, by choosing properly the molecular energy parameter J as a function 
of A; in this way the critical temperature of the system with A = and A 7^ is the 
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same. We consider then the model (1.1) with Jr replacing J, and we will choose 
Jr = J -\- 0{\) so that the critical temperature is in correspondence tc = — 1. 
We define 

< 0(x)0(y) >A= ^ Yl 0(x)0(y)e-^-('^^'^'-''^) (1.10) 



where Z\ — ^^(i) ^(a)^.^^ e ^^^^^ ^''^^ is the partition function. 



f^{J) = lim|A|^oo -pj log^A is the free energy, -§jif^{J) is the specific heat, and 
the correlation function of the observable 0(x) is 

< 0(x)0(y) >A,T=< 0(x)0(y) >a - < 0(x) >a< 0(y) >a • (1.11) 

In the Ising model close to the critical temperature -§jif^{J) ^ Clogjt — tc\~'^ 

and for large distances | < 0(x)0(y) >a,t | < |x-y|^' — ~"' '^^^^ A,C suitable 
constants. 

We will prove the following theorem. 

Theorem 1. Consider the hamiltonian (1.1), with Jr replacing J, and assume 

that \a -\- h\ > and (1.9). There exist £ > 0,a3 > 0,a > constants and a 

_ 1 

function i/(A) = tanh — tanh J such that, for |A| < e and a > \t — tc\ > e , 
the energy-energy correlation verifies 

hm < ct4«)4?V(")4?) >A,T=^^(")'"(x,y) + ^(")'^(x,y) + ^(")'^(x,y) (1.12) 

|A|^oo 

where x, x' and y, y' are nearest neighbor pairs respectively and, for any integer 
positive N 

where Cn,t are positive constants, A = \t — td^'^'^^ , and ^{X) = 0(A), 

r7i(A) = -aiA + 0(A2) ?72(A) = aaA + 0(A2) r73(A) = aiA + 0(A2) (1.14) 
with ai > 0, a2 > constants, i/ 1 < |x| < 
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and |n^"^''^(x, y)| < |^_y|2+T ? where Zq > is a constant and /^(x, y) are smooth 
bounded 0(A) functions. Finally the specific heat verifies 

Ci|-[l-|A|2^i] + l[l-|A|2''3]| < \-^f\j)\ < C2\-[l-\A\^'''] + -[l-\An 

(1.16) 

where Ci , C2 are positive constants. 

1.6 Remarks. The above theorem describes the behaviour of the energy-energy 
correlation and the specific heat near the critical temperature of Ising models in 
two dimensions weakly coupled by a four spin interaction. One can distinguish 
two different regimes in the asymptotic behaviour of the energy- energy correlation 
function, discriminated by an intrinsic correlation length ^ of order \t — tc\~^~'^'^ 
with r]2 = 0{\). If 1 << |x — y| << ^, the bounds for the correlation function 
is power-like (with a A dependent exponent) while if ^ << |x|, there is a faster 
than any power decay with rate of order In the first region we can obtain 
the exact large distance asymptotic behaviour of the energy-energy correlation 
function, see (1.15), while in the second region only an upper bound is obtained. 
The logarithmic behaviour of the specific heat in the Ising model getting closer 
and closer to the critical temperature is changed by the four spin interaction in a 
power law, at least up to temperatures very close to tc for small four spin coupling. 
Abound like Ci log(|t - tj"^) < \-^f^{J)\ < C2 log(|t - 1^"^), for some positive 
A-independent constants Ci and C2, which would be true if the model would be 

in the same universality class as the Ising model, cannot be true by (1-16) at least 

_ 1 

for \t — tc\ > e ^3>^. 

In the particular case a = c — our model reduces to the Eight-vertex model, 
which is exactly soluble [Ba]; our results agree with the informations obtained 
by the exact solution, in which non universal critical behaviour is found and it is 
believed that the specific heat diverges with a power law [W] . In the case a = b = 
our model reduces to a model of two non interacting Ising models with nearest 
neighbor and four spin interaction within each copy, which was studied in [PS] 
(which is indeed the first paper in which fermionic RG methods were applied to 
classical Ising-like models, and it is a major source of inspiration for the present 
work). In this case Cilog(|t - tc\~^) < |^/^(J)| < Oslogdt - tc\~^) (up to 
t = tc) so that universality indeed holds 

1.7 Sketch of the proof. The proof of the theorem is based on Grassmann vari- 
ables combined with renormalization group techniques. In §2 we briefly recall the 
well known representation of the bidimensional Ising model in terms of Grassmann 
variables, mainly due to [K], [H], [MW], [S] and recently rederived in a more cohe- 
sive way in [PS]. The partition function and the correlation can be written in terms 
of Grassmann integrals with a quadratic action which can be explicitly computed 
in terms of Pfaffians. The exact solvability of the Ising model is hence related to 
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the fact that it can be expressed in terms of free fermions, as it was first noted in 
[SML]. In the Grassmann representation there are four independent Grassmann 
variables associated to each point of the lattice, two with a large 0(1) mass and 
the other two with a continuously vanishing mass at the critical temperature. In 
§2 we show that two Ising models weakly coupled by four spin interactions, with 
hamiltonian (1.1), can be also written in terms of Grassmannn integrals, with the 
difference that there are now eight independent Grassmannn variables (four for 
each Ising model) and that the formal action contains now terms which are quar- 
tic in the fields, corresponding to a short range interaction among fermions. By 
performing a suitable linear transformation in the Grassmann variables one ob- 
tains that in the new variables the quadratic part of the action strongly resembles 
the action of two massive Dirac fermions in d — 1 + 1 dimensions (it would coincide 
with that if the continuum limit would be taken); again one Dirac fermion has a 
large 0(1) mass and the other is vanishing at tc, and we call them heavy and 
light fermions. We are essentially exploiting in this transformation the relation 
PfA = detA. 

In §3 we integrate out the heavy fields in the Grassmann integral for the partition 
function, so obtaining a Grassmann integral whose formal action contains mono- 
mials of every degree in the light Dirac fermions. In order to do this we use the 
representation in [Le] of fermionic truncated expectations and Gram-Hadamard 
inequality. 

In §4 we apply renormalization group methods to integrate the light fermions, 
which in a sense are the critical modes. We will use a suitable modification of the 
multiscale expansion used in [BM] (see also [BG] or [GM] for a general introduction 
to the formalism) to study the correlation functions of the Heisemberg- Ising XYZ 
chain; the close relationship between Ising models with quartic coupling and the 
XYZ times has been pointed out many times in the literature, see for instance 
[Ba]. A power counting analysis says that the terms bilinear in the Grassmann 
variables are relevant in a RG sense, while the quartic terms, or the bilinear with 
an extra derivative are marginal. One can understand here why universality is 
still present if one considers decoupled Ising models with nearest neighbor and 
four spin interaction within each copy, and why on the contrary it is lost if the 
Ising models are coupled. In the first case (which is the one treated in [SP]) 
one can easily check that the local part of the quartic terms is vanishing, so the 
effective interaction is indeed irrelevant in the RG sense. On the other hand in the 
second case the quartic interactions is truly marginal and this produces a line of 
fixed points for the RG transformation (instead of the gaussian fixed point as in 
the previous case) continuously depending on the coupling A. We decompose the 
Grassmann integration P{d'^) as a product of independent Grassmann integrations 
P{dilj^^^), with covariancc with non vanishing support only for momenta with 
modulus between 7^^"^ and 7'^''"^, with 7 > 1 and = 0, —1, —2, —3, .... We 
integrate each P{dip^'^^) iteratively starting from P{dip^^^) obtaining a sequence 
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of effective potentials V^'*^ describing the theory at momentum scale 7^*; at each 
step new contributions to the mass and the wave function renormalization are 
obtained which are included in the fermionic integration; hence P{dilj^^^) has a 
covariance with mass rrih and wave function renormalization with a non trivial 
dependence on h, i.e. ruh — \t — tc\j^^^ and Zh ~ ^'^4/1^ with 774 = O(A^). The 
iteration stops as soon as the mass m/j becomes of order 7^*; we will call h* the last 
scale to be integrated (of course h* — > —00 at the critical temperature t = tc). The 
iterative procedure allows to write the effective potential V'^'*-* as sum of monomials 
in the Grassmann variables, with coefficients which are (convergent) perturbative 
expansions in terms of a few running coupling constant, Xh, the effective coupling of 
the interaction between fcrmions, Uh, which takes into account the renormalization 
of the value of the critical temperature, dh related to the renormalization of the 
fermionic velocity and other couplings Ci,ui,h, Ci,oj,h which are the coefficients of 
quadratic terms in the Grassmann fields with a derivated field. Despite there 
are many similarities between our Grassmann integrals and the one describing 
relativistic fermions, there is a crucial difference; the interaction term in the action 
is not invariant under gauge transformation, hence terms which were absent in the 
free action can be generated in the RG iterations. One can check that indeed 
peculiar symmetries of the model (1.2) ensure that there is only one relevant term 
quadratic in the fermions (if there were more one gets in troubles, as there is only 
one free parameter in the hamiltonian). On the other hand we cannot exclude the 
generation of marginal quadratic terms which were absent in the free action; they 
are the terms Ci,Lo,h,Ci.Lo,h- At the end the result of this iterative integration is 
an expansion for the partition function in terms of the running coupling constants 
which is proved in §4 to be convergent provided that the running coupling constants 
are small for any h. 

In §5 we prove indeed that it is possible to choose the counterterm i/ as a function 

of A so that the running coupling constants are indeed small; the condition on the 

_ 1 

temperature |t — ^ e is used to control the flow of the running coupling 
constants, and some cancellations at the lower order in our expansion are also used 
to be as close as possible to the critical temperature. 

Finally in §6 we define an expansion for the correlation functions and the specific 
heat; it is similar to the one for the partition function, with the main difference 
that one has to introduce new fields associated to the external fields. There are 
additional marginal terms in RG expansion to which other renormalization con- 
stants, with a non trivial behaviour in h, are associated, i.e. zj^^ ~ 7^^'* and 
zj^) ~ ^'73/1 _ By guch expansion the statements in the theorem are derived. 

Among interesting open problems there is the reaching of the critical tempera- 
ture; surely one has to exploit suitable cancellations at every order of the expansion 
for Xh,5h, as in the theory of d = 1 interacting Fermi systems, see for instance 
[BGPS] and [BM], but at the moment the main difficulty is in the fiow of the cou- 
plings Ci,uj,h and Ci,uj,h- Another interesting problem would be the analysis of the 
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spin-spin correlation function whose expression in terms of Grassmann variables 
is unfortunately quite complicated and not easy to manage. A similar problem 
appears in considering two Ising models coupled by a weak interaction bilinear 
in the spins; in such a case the interaction in terms of Grassmann variables has 
a quite complex expression and it seems difficult to study. On the other hand, 
for large coupling, such bilinear interaction should be irrelevant and universality 
should hold. Finally it should be interesting to study the case of two coupled Ising 
models at different temperatures or the case of four coupled Ising models; in this 
last case interacting spinning d = 1 fermions appear in the fermionic description, 
which are known to have a behaviour quite different from the spinless one (like in 
the d — 1 Hubbard model). 

2. Fermionic representation 



2.1 Grassmann integrals. If A = the hamiltonian (1.1) is given by the sum of 
two independent Ising model hamiltonians, and the partition function is given by 



Z^j^^zP where 



Zj'^)= J2 e-^^(-^"'). (2.1) 



It is well known that such partition function can be written in terms of Grassmann 
integrals, so we recall first their definition. Grassmann variables r]a, a = 1,2, .., 2n, 
,n even, are anticommuting variables satisfying 

{Va,Va'}^0 77^=0 (2.2) 

The Grassmann integration J dr]a is a linear operation defined as 

dr]a=0 / dr]ar]a = 1 (2.3) 



It holds that 



a 

where A is an even antisymmetric 2n-matrix and PfA denotes the Pfaffian. It 
holds 

J <ir/i...dr/2n cxp i ^ VaAa,pVf3 = (2.5) 

a, 13 
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where the sum is over all the permutations. We can consider another set of Grass- 
mann variables rj^ , a = 1, ..,2n, and 

[Udvcll dv+ exp[^ VaBa,i3V^] = DetB (2.6) 

a a a, 13 

The well known relation (Pf^)^ = det^ can be quite easily deduced by the above 
Grassmann integrals; it can be written as 

/ Y[dr^^dr^U''"''-'^'^t = j J]d77ii)e^''«'^-^<' j J] c^^i'^e^"" (2.7) 

which can be proved by the change of variables 

rit = ^(ryi^) + = ^(^i'^ - iri^a^) (2-8) 

in J drjadrj^e^"^"''^^!^ . Then drjadrj^ = idrja^drja^ and 

VaAa,pvt = Iv^a^A^W^^ + Iv^a^ A^,Pvf (2-9) 

as 

V^a^A^,f,rjf - r/i^)^a,/3r/W = rjl^^A^,fsvf - V^^^ Af,,^rjS^ = (2.10) 

2.2 The fermionic representation of the Ising model. It is well known that 
the partition function can be written as a Grassmann integral. It is proved in 
[K],[H],[MW],[S] that 

Zj") = (coshJ^)^2^^ f Yl dHi''UH^"Uv^"Uv^''\-e^+'++e'^+'-+e^-'++e'^-' 

(2.11) 

where B is the total number of bonds and S is the total number of sites, 

Y: tanhJ.[ifH„ifiJ,,,„ + V;(^^<t+i] + 



Elffio:) TT{a) I T/(o!) TA(a) i T/(o!) tt{oi) _i_ t/(o!) TT{a) < zria) yia) lyia) tt{ 
I X,Xo X,Xo ' X,Xo X,Xo ' ' X,Xo X,Xo ' ' X,Xo X,Xo ' X,Xo X,Xo ' ' X,Xo X 



X,Xoi 

x€Am 

(2.12) 

and Fi") are Grassmann variables such that 

x,xo+M '^■^■^x,xo ^■^x+M,xo x,xo 

^x,xo+M — ^^x,xo ^x+M,xo ~ ^ ^x,xo y^.±0) 
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and identical relations hold for the variables V^°'^ , V^°^^ . The Grassmann integra- 
tion 

/ ria; du'^^ dll'^^ is defined as the linear functional on the Grassmannian algebra, 
such that, given a monomial Q{H^'^\ H^"^^) in the variables H^\h^\ x e Am, its 
value is 0, except in the case up to a permutation 

of the variables. In this case the value of the functional is determined, by using 
the anticommuting properties of the variables, by the condition 

W H^^H^^ = 1 . (2.14) 

xeAm 

In a similar way is defined the Grassmann integration for V^°'^ , V^°'^ , just ex- 
changing H,H with V,V. In order to prove (2.11) the starting point is the high 
temperature expansion 

Z}"^ = (cosh Jr)^2^ ^(tanh J^)' (2.15) 

where the sum is over all the closed polygons which may have points but not sides 
in common, see Fig. 4. One then replaces the Ising lattice with another lattice, 
called second lattice, in which each site is replaced by four surroundings sites or 
terminals, see Fig. 5. If we associate with AT^., with S'x, -ffx with and 
with Wx, see Fig. 5, it is easy to verify, by the rules of Grassmann integrations, 
that / e^-.-' can be represented as a sum over polygon configurations over the 
second lattice. Such polygon configurations are defined such that each terminal 
must be coupled to another terminal, Ey^ may be coupled to W^^i^x^; ^x„,x may 
be coupled Sx,xo+i^] may be coupled to iVx, Wx, .S'x; A^x may be coupled to 
W^x or S'x; VFx may be coupled to S-x., and no other types of couplings are allowed. 
An example of polygon configuration in the second lattice is in Fig. 5. By a local 
identification of Ising configurations and configurations over the second lattice 
(see Fig. 6), there is a correspondence between polygon configurations on the Ising 
lattice and polygon configurations on the second lattice. The correspondence is 
not one to one, as the configuration with zero lines in the Ising model corresponds 
to three configurations in the second lattice. However the signs add so that the 
sum over polygon configurations in the Ising lattice in (2.15) is exactly equal to 
the sum over polygon configurations in (2.11); the proof is based on a subtle 
combinatorial analysis first done in [K]. In order to ensure periodic boundary 
conditions, one has to sum over four terms as in (2.11), as explained in [MW]. We 

write e = e e where S). ^! is given by (2.12) with J replacing and 



.xeA» 
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If J is not constant but it depends on the bounds one obtains a similar formula 

e,e' 



in which S^f^l'^ is given by 



e,e' ~ 2-^^'^^^^^'^1;x,xo;x+1,xq^x,xo^x+1,xq ^ ^'^^^^^ 2;x,xo;x,xq + 1^ x,xo^ x,xo+l\^ 

(2.17) 

ErrT(a) rr(a) , T/(a) T/(a) , y{a) TT{a) , T/(a) o-(a) , o-(a) T/(a) , yia) TT{a) 1 
1 X,Xo X,Xo X,Xo X,Xo ' X,Xo X,Xo ' X,Xo X,Xo X,Xo X,Xo ' X,Xo X,Xol 

X 

and the factor (cosh J)^ is replaced by ]^(^cosh Jf,, where the product is over all 

the possible nearest neighbor bounds. We will call Zj"''(Jx,x') the Ising model 
partition function with non constant J. 

2.3 The fermionic representation of coupled Ising models. The partition function 
of (1.1) is 



Z2i= J2 E e-^^(-^ ^)e-^^('^^ ')e-^('^^ ^) (2.18) 
Let us consider first the essentially equivalent expression 



a(l)=±l ,7(2)=±1 



X X 
X X 

UUi^+^<^^i%^it,xo<^i%<^il^^^^ 
ax ax 

Noting that 

(a) (a) _ )l r ^ (2 20) 

(^x,Xo'^x+l,Xo^ - (a) ^/ '''''' ^'{^^"i,}={-^'"'} 

'^'^l;x,xo;x+l,xo 

and from (2.17) this derivative gives an extra factor tanh J*^")+sech^ J*^'^)^a%^oiy^^\ .^.^ 
in (2.11). In the same way 

"^x.xo^x.xo+l^ - ^^(a) x,x';'l{jM (^Z.Zi; 

'^'^2;x,xo;x,xo+l 
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and this derivative gives a factor tanh J^'*) + sech^ J^'^^vi^ao^i'^o+i- We can write 
than, if 5+^+ = 1 and = (^-,+ = (^-,- = 2 

Z2z= Yl h^f^'''-'''' E {-l)'^''''^''''Zl'j ''''''' ''''''''''' (2.22) 

e(i),£'(i) e(2),£'(2) 

where 

^e(^),s'(^),£(^),s'(^) ^ (cosh J)2^22^i (2.23) 

a=l X 

+ Aa(tanh J + sechV^W„ifWi,,J(tanh J + sechV^gl^if^i^.J] 

X 

+ Aa(tanh J + sechVv;^ V;« +i)(tanh J + sech' JV^%v}%^,)] 

X 

+ A6(tanh J + sechV^W^i^i^^^i,, J(tanh J + sechV^ifi^ vjS^+i)] 

X 

+ A6(tanh J + sechVV.W vi;i^+i)(tanh J + sechVi72i,.o+i<io+i)] 

X 

n nil + Ac(tanh J + sechV^(^iiy(Ji,,J(tanh J + sech' JV^^J^V^'^J^^,)] 

a X 

n nil + ^^(t^^h + sechVi^^-iiy^Ji,, J(tanh J + sech' JV^l{,^_,V^llj] 

a X 

The above expression can be rewritten as 

^.(^),£'(^),s(^),s'(^) = (cosh J)2^22^i /" [H n^^x"^rf^4"^rf^J"^^^^i"^e^^""^-'^"^]e^ 

a=l X 

(2.24) 

with 

V = Va + Vb + Vc (2.25) 
and, if fi = log{l + X[i] tanh^ J) and [i] = a,b,c 

X 

+ E[/« + ^«[^i2o<'io+l + ^ifio^l+l] + Aat^i'VW^+it^i3„Ki2„ + i] (2.26) 

X 

X 
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(2) 1 

X,Xo + li 



X 

X a 

I Y^Y^Ff _|_\ fT/C") T/(") 4_H-(a) l4_\ t7('*) fj{cc)TT{oi) 

^ ^Ljc-rAcL''a;+l,a;o-l '^a;+l,a;o"'"-"a;,a;o-"a;+l,xoJ"'"^c •'x+l.xo-l ''^a;+l,a;o-"x 
X a 

It is easy in fact to verify that 

(2.27) 

hence the equality between (2.23) and (2.24) holds with the identification 

Ai(l + tanh^ J) = A[z]sech^ Jtanh J 

(1 + tanh^ J) (A, + (A,)^) = A[i]sechV (2.28) 

An expression identical to (2. 22), (2. 23) holds for Z21, the only difi'erence being 
that the relation with respect Aj, Aj is slightly more complicated than (2.28), but 
for small A again 

\i = [z]A(tanh JsechV + 0(A)) Aj = [z]A(sech^J + 0(A)) (2.29) 
2.4 Massive and massless fermions. We define 

PsS) £'(«) (^^''' ^^^"^) = n f^^^4''^f^^4'*^f^^i"^'^^i''^e^^^"^'^''"^ (2.30) 

X 

It is convenient to perform the following change of variables [ID] , [SP] 

Fi") + = + xi"^ K^"^ - ^^i"^ = + xi"^ (2.31) 
It holds that 

Plt^.^dH^^dV^''^) = ^'iS,,'(«)(#'^)^^S,,'(«)(^x'^)e^^^'"''^'"'^ (2.32) 
where, if t — tanh J 



P«,,(#(«)) = [ n di;t^di;^^^]eMl E V'4'^)(5i-^5o)V'4'^)+V54'^n^i+*5o)V54'^^) 



-4 
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+ J2 ^(V2-l-^)4"V4"^} (2.33) 

where 

dii^t^ = 4t,.o - ^^"^ ^oV'i"^ = V'il+i - V'^"^ (2.34) 

Moreover 

^ E [-^xi"H5ixi"^+5oxL"^)+^xi"H5ixi"^+aoxi"^)]- E ^(V2+l+^)xi")xir^} 

x^Am xSAm 

(2.35) 

and finally 



X 



X^^\^^^|;i-^ + i^o^|Ji^^) - X^'^H^iV^^"^ - ido^l^i^^) + #(5iX^"^ - ^oX^"^) (2-36) 
+i^t\-diX^"^ - dox^"^) + X^^H-SiV'^"^ - doi^i"^) + X^^H-dii^t^ + doi^t^)} 

We find convenient to rewrite the Grassmann variables in momentum space. We 
call De,e' the set of k such that 

_ 27rni {e - l)7r _ 27rno (^^ - l)7r 

and -[M/2] < no < [(M - l)/2], -[M/2] < m < [(M - l)/2], no,ni e Z. We 
can write if k = {ko, k) 

Hence 

n^k#) = [ n ^^"^^V'i'^^] expl^^ E [V'i'^VS(^sin/c + sin/co) 

+4"^V'-k(^sinA; - sin/co) + i2m^(k)V;i"V-k]] (2.39) 
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where m^(k) = f (^2 - 1 - t) + (2cos/co + 2 cos /c - 4) = 2^^^ + 0(k2). In 
deriving the above formula we have used that 

^V^^^aoV'^"^ = J2 (e^'^" - l)V'i"Vik = I] isin/coV'i^V-k (2-40) 



as 

5] (cos/co-l)V'i"VS = - E (cos/co-l)V'iVi"^ 



= - J] (cosA;o-l)V'i"Vi'£ = (2.41) 
2.5 Majorana and Dirac fermions. We write ■0x"'' in (2.38) as 

where = (^^^J^, ^^w^), and •i/'J"^ = '^k"ip, moreover 



pj'^] (d^) = p(«) (dV'Oe^'^"^'^' (2.43) 



s(«) 



e,£ 

where 



k'eD_,_ k'eD_,_ 

(2.44) 

V^-k' (^ sin /c' - sin k'o) + i2m^ (k^V^k^" Vl^^' ]] 

and (5^ e £/ is defined by (2.43) and is formally vanishing in the limit M — > oo. 
Proceeding in the same way for P{dx) and Q we find 

/ [fl ntle'(«)(^V'")nt) (^X'^)]^^^^''^^^''^^'^^ = (2-45) 



a=l 



a=l 



^Qe{^l>',x')^Q{^',x')^V{^',x') 



where 

2 



4(V'',X') = E[^£s(«),e'(«) (2-46); 
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of course Q-,-,-,-(V'')X') = 0- If 

= {k' e : A;' > 0} (2.47) 

we can write P^"-\dip') in the following way 

k'eD+ _ k'eD+ _ 

(2.48) 

&^l = {i^'i^\i^S^.^'i^\^'S') (2.49) 



where 



and ^(k') is 



/ sin A;o + z sin A;' 

- sin A;q — z sin A;' — ■imi/,(k') 

zm^(k') — sin A;q + z sin A;' 

\ sin A;o — i sin k' 

(2.50) 

Hence if we perform the change of variables 

Vi> - ^{^i^^ + i^S^) V'lV = ^(V'k^'^ - i^S^) (2.51) 

we find by (2.7) 

pW(dV')^^^n#') = (2-52) 



n n #^,.#k',c.#^k',.#:k',c.exp[^ J] f^'A(k')^lt)] 



where 



^k = (V'k,i' ^-k,i' ^k,-P ^-k,-i) (2-53) 

f^'^ = «i,V'^k,i'V'J-i'V'Vi) 

We can make another change of variables i^^^ — > V'^k uj ^^^^ 

P{d'^) = (2.54) 



n n ^V'^k,.#k,.#j.#:k,.exp[^ 5^ ijA(k)ei 



c(+)l 



'2M2 

where 



^k = (V'k,i'V'-k,i'V'k,-i'V'-k,-i) 
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ii''^ = (V'lk,i,<i,V'lk,-i,V'J-i) (2.55) 
and in terms of the original variables 

^1,^' = -^ii^'^'^ + ^^'S^) ^it-k' = ^(^'i^^ - (2-56) 

In the following we will call k' simply k. If we write 

^^ = W^ ^^''^V'J (2.57) 

keD_,_ 

we can rewrite (2.56) in the coordinate space as 

V'l^x = ^(^?^ + = ^i^'^'^ - ^^'^'^) (2-58) 

V'li.x = ^(V5?) + #)) V'ii,x = ^(4^^) - 

or in terms of the original variables 

= -^(e^^e(^).'^^)^V4^) - ie^Pe(^).'(^)-V'4^)) (2.59) 

V^l,,, = -L(e^Pe(l),e'(l)-Vi4^) _^e^Pe(^),e'(^)-Vi4')) 

If k e then -k e i:*" with D+ \\D~ = D_ _ so that 

P(#) = (2.60) 

n n ^V't.#ir,.exp[^ ?TA(k)et+)] 

keD_,_u;=±l keD_,_ 

if-L-^ _ i sin /c + sin /cq im^(k) \ 
^ y — zm^(k) z sin — sin A;o / 

= (V'k.i, V'k,-!) fJ'^ = (V'J 1, V'J-i) (2.61) 
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Remark. In the physical language, the change of variables (2.59) means that one 
is describing the system in terms of Dirac fermions instead in terms of Majorana 
fermions. 

A similar computation can be done for P{dx)', proceeding exactly as above we 
find 

p(i)(^^(i))p(2)(^^(2)) ^ p(^^) ^2.62) 



where 



P(dx)= n n <-^Xk,.exp[^ Yl CB(k)4+^] (2.63) 



ker> u!=±i keD_ 



\ —im-^{k) z sin A; — sin A;o / 

C = (Xk,i' ^7+'^ = (x£i, _i) (2.64) 

Note that tm^(k) = 2(v^ + 1 + 1) + 2t{2 — cos /cq — cos ki), and the mass of the x 
fields never vanishes. It holds that 

/ P{dx) xx,c.x+.' = g^^M - y) . (2-65) 



where 



,(x) 



^.vl'^-y) = , (2.66) 



k 



and -B-i(k) is the inverse of the S(k) defined in (2.63). 
If we set 

A(y) = det S(k) = - sin^ ko - sin^ k - [m^(k)]2 , (2.67) 

then 

B-'^(\A= ^ A-sinfcp +isinfc -im^(k) \ gg) 
^(k) V im^(k) sin /cq + « sin /c J ' \ • ) 

Similar formulas hold for g'^'^\-x. — y) and, if i = V'jX) the following bounds holds, 
for any A'^ > 1 

ie(x - y)l < i^|d„V-y)|i + Kd!(.-,)|. 



l''"-"<''-''>l-l + Kd„(x-y)IN (2-™) 



where 

dM(x - y) = (- - «^^( — M — ^ ^^-^^^ 
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The same transformations are done for the term e which can be written as 
the product of three terms similar to (2.33), (2.35) and (2.36) with t replaced by 
2^ and v/2 - 1 in (2.33) and v/2 + 1 in (2.35) replaced by 0; then (2. 33), (2. 35) can 
be written in terms of Dirac fermions as (2.61) and (2.64) with t replaced by u 
and \/2 — 1 or \/2 + 1 replaced by 0. 

We have then written the partition function as Grassmann integral over the 
fields, which are massive with a mass 0(1), and the V'-fields with a small mass 
0{t — tc) close to the critical temperature; we integrate out first the heavy fermions 
X to get an effective theory in terms of the light fermions only. 



3. Integration of heavy fermions 

3.1 Local interactions. By the change of variables in the preceding section we 
can write (2.23) as 

where jV is a normalization constant, P{ip) and P{dx) are given by (2.60), (2.63), 
V, Q, Qs are obtained respectively from (2.36) and (2.46) by the change of variables 
(2.59). We write 

V(V, X) = V2(V) + V4(V') + V^(V', x) (3.2) 

where V2{ip) is a sum of monomials bilinear in the ip, ip'^ variables, V4{'ip) is a sum 
of monomials quartic in the and V^{il), x) is sum of monomials bilinear or 

quartic each one containing at least a Xuj field. It holds that, from (2.31), 

H^^^H^^^ = '^i^t^i^t^ + n-^v^"^ = '^i^t^i^t^ + Rt^ (3.3) 

where -R^""* is sum of monomials bilinear in the fields %, x, ■0, '0 and containing at 
least one field x o^" X- Moreover, from (2.58) 

i^'^H^^^ = + + V'ltxV'-i.x + V'JxV'li,x] (3.4) 

V'^'^V'^'^ = ^[-V'i,xV'-i,x - V'lVii.x + V'iV-i,x + V'i,xV'ii,x] 

hence 

^(1)^(1)^(2)^(2) ^ e-^^(Pe(^),e(l)+Pe(l),e(l))'^^+^^+^_^^-^^^^-^ (3.5) 
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and from (2.26) 

V,{^) = -2e-^^(P^(^).^(^)+P^w,^w)-(A„ + A,) J] ViVii,xV'Ii,xV'i-x + (V) 

(3.6) 

and is a sum of quartic monomials with coupling 0(A) in which at least a di(^ 
field. In the same way it is easy to check that 

"^2 = 1^ + /i(A)][ia; + cos/co + cos /c - 2]V'k,a;V'k 



k. 



UJ 



[u + /2(A)] [a; sin ko + i sin k]i;k,u>i^+J (3.7) 
where |/i|, I/2I < C|A|. 

3.2 Integration of the heavy fields We now integrate the heavy x fields 

j P(#)e^'^'''+^'''^^) = j P(#) j P(dx)e^(^'^)e'5^(x,^)eV(^,x) (3.3) 

where jV^-*^) is a constant. Let us consider a set of coordinates xi, X2n which are 
not all different one to the other and we will denote by V* ^ the sum over 
all the distinct < 2n variables. We prove the following result. 

Theorem 2. There exists an e such that, for |A|, jivj < s 

* 

^^'^ = E E E ^n(xi,..,x2.)a«^V'^;,,^...a-"V'S-,,^„ (3.9) 

n>l {s,a,ui} xi,..,X2„ 

and, for n > 2 

* 

J2 |t^n(xi,...x„)|<M2C-|Ar/2 (3.10) 

Xi,...X2„ 

For n = 2 

* 

^ w^2(xi,..,x2n)a«^V^;,^,...a«-v^-,,,„ = 

Xi,..,X4 

* 

E E w^2(xi,..,x4)a-^V'Sl,c..-5"^"V'^r„,u,2„ (3-11) 

Xi ,..,X4 {g,a,Ci;} 

Exi,..,x4 \W2{^i, X4)| < M^C\X\ and W2{0, 0) = ©(A^) ; /orn = 1 



^ J2 V(^i(xi,X2)a'^^Vxi''^'^'*'V': 

xi,X2 {e,a,(t;} 



£2 

X2,a;2 
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M^)+Y^ E ^l(^l'^2)^"^^xl,c..^"^V'SL,c..„ (3-12) 

with Exi,x2 l^i(^i'^2)| < M^C\X\; Vi and V4 are given by (3.6) and (3.7). 
Finally 

P{d^)= n n ^V'J.#k,.exp[-^ V'J.i.,.'(k)Vkx] (3-13) 



with 



i(k) 



Z{isink + sin/co) + /i,i(k) -i{t - tc)Co + /i,2(k) 
Co + /o,o(k)V *(^-*c)Co + /2,i(k) Z(zsinA:-sinfco) + /2,2(k) 



TO! 



t/i Co = (t + 1 + v^)2 and Z = |[(2t + 2x/2t) + (2^2 + 3 + t^)] and /o,o(k), 



/ij(k), i,j = 1,2 analytic 0(k ) functions. 

Remark 1. The fact that M^2(0, ...,0) = O(A^) can be checked by an exphcit 
computations of all the contributions with coupling 0(A) to W2, noting that they 
can be only obtained contracting a terms quartic in the x fields with one of the 
addend of (2.36); each of such terms carries a derivative in the coordinate space, 
hence the Fourier transform of such terms is vanishing at zero momentum. 

Remark 2. The fermionic integration P{d'i(j) reduces, neglecting the functions 
/o,o(k), /ij(k) which are O(k^), to the integration of a system of Dirac fermions 
on a lattice with mass 0{t — tc). Hence, apart from the functions /o,o(k), fij(k), 
it coincides with the free action of the Thirring model. Note however that the 
interaction in the Thirring model is quartic and invariant under the transformation 
V'x.o; ~^ e*^'*'^Vx,a' ^fi® massless case; both this properties are not true for (3.9). 

Proof. We start from the definition of truncated expectation: 

^?(^; ^) = ^ log J Pidxy^""^ |a=o (3.14) 



so that, calling 
we obtain 



V(x, ^) = Q(x, ^) + Qe{x. ^) + V(x, ^) (3.15) 



n=0 

We label each one of the monomials (whose number will be called Co) in V by an 
index Vi, so that each monomial can be written as 
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where x„. is the total set of coordinates associated to Vi and Py^ and P„. are set 
of indices labehng the x or V^-fields. We can write 

VW(V')= Yl (3.18) 

v<"(p.,) = El n s°"'CfA-(/)i4!>".) . (3-19) 

oo ^ n 



'(0) 

y-^voj — 

n=l vi,..,Vn i=l 



4!;(^-o) = E;^ E s^[xiPvj,---.xiPvj]ii^iM, (3.20) 



where Z)a,i,...,;„ < Q > -Pfo = Ui -Pi^i and = \J- x„.. We use now the well known 
expression for (see for instance [Le]) 

£j(x(Pi),...,x(P.)) =En^^,^+(^'-y') / ^^^T(t)detG^(t) (3.21) 

where: 

-)P is a set of indices, and 

X(P) = n d-'''x^S)Mf) (3.22) 

-)T is a set of lines forming an anchored tree between the cluster of poins Pi, ..,Pg 
i. e. T is a set of lines which becomes a tree if one identifies all the points in the 
same clusters. 

-)t = {ti^i' e [0, 1], 1 < < s}, dPrit) is a probability measure with support 
on a set of t such that tj^j/ = Uj • Uj/ for some family of vectors e of unit 
norm. 

is a {n — s + 1) x (n — s + 1) matrix, whose elements are given by 
^ hi'9u^-,u^+i^ij - yi'j') with iffjj^j,) not belonging to T. 
If s = 1 the sum over T is empty, but we can still use the above equation by 
iterpreting the r.h.s. as 1 if Pi is empty, and detG{Pi) otherwise. 
We bound the determinant using the well known Gram-Hadamard inequality, 
stating that, if M is a square matrix with elements Mij of the form Mij =< 
Ai, Bj >, where Ai, Bj are vectors in a Hilbert space with scalar product < •, • >, 
then 

|detM| < • . (3.23) 

i 

where || • || is the norm induced by the scalar product. 

Let Ti. — M.^ (E) H-o, where TYq is the Hilbert space of complex four dimensional 
vectors F(k) = (Fi(k), . . . , F4(k)), Fj(k) being a function on the set I'-,-, with 
scalar product 

<F,G>=J2^T. ^*(k)G'.(k) . (3.24) 

i=l k 
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and it is easy to veriy that 

(3.25) 

where G M^, i — 1, . . . , s, are the vectors such that ti^i' = Uj • Uj/, and 

4 - ik'x 1 f (-sin/co +isin/c,0, -zm^(k),0), if w = +1, 

"'"^ ^ " ■ I (0' ^^x(k), 0, m^(k)), if c = -1, 

R _ ,ik'y 1 1(1,1,0,0), ifa; = +l, 

y=Z(k) \(0,0,l,(sin/co + isin/c)/m;,(k)), ifa; = -l. 

(3.26) 

Hence from (3.23) we immediately find 

|G5,,,,.|<Cr (3.27) 

Note that Ci is an 0(1) constant; this foUows from the fact that the Grassmann 
variable x h^s an 0(1) mass. Finally we get 

oo ^ n 

Ei4li('^-)i^E;^ E E ^rE[ni^x('^'-y')i]ni^^('^-)i 

x„o n=l ■ x^^ ,...,x^„ T leT i—l 

(3.28) 

where we have used that J dPrit) — 1. Regarding the sum over T, it is empty if 
n = 1. If n > 1 the number of anchored trees with di lines branching from the 
vertex Vi can be bounded, by using Caley's formula, by 

\PvA''...\PvJ''-, (3.29) 



(di-l)!...(d,„-l)!' 



hence the number of addenda in is bounded by n!C^. Finally T and the 
[J^ x^. form a tree connecting all points, so that using that the propagator is 
massive and that the interactions are short ranged X^TfllzeT Idxi^i ~ 

yi)\] nr=i < C3 |A|'^M^, where n is the number of coupling 0{X). 

Let us consider the case \Pvo\ > 4. Note that if to Vi are associated only terms 
from V(V',x), th^^^ ^ = The same bound holds for M large enough, if there 
are Vi associated with Q^; in fact such terms are vanishing as M — > 00 (one has 
in the bounds an extra for any of such Vi, for dimensional reasons). Let us 
consider now the case in which there are end-points associated to Q{i(^,x)j which 
have 0(1) coupling; there are at most \Pvo\ end-points associated with Q{ip,x)- 
In fact in (5(V', x) there are only terms of the form i/jx, so at most the number of 
them is equal to the number of fields. If we call nx < "h the number of vertices 
quartic in the fields it is clear that n\ > max{l, \Pyg\/2 — 1}; hence 

00 

y (x„J| < M^y C^+l^-ol|A|t|Ar^^{^'l^^ol/4-i/2} (3.30) 
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and (3.10) holds for \Pvo\ > 4 (in the r.h.s. of (3.30) the sum over the number 
of I'i to which are associated quadratic monomials with coupling 0(A) is already 
done) . 

Consider now the case \Pvo\ = 2; in this case there are terms A independent, 
obtained when to all the Vi are associated with elements of Q{'il', x). It is convenient 
to include all such terms in the free measure, as they cannot be considered as 
perturbations (they are not 0(A)). Instead of computing all such terms, we can 
proceed in a more rapid way by noting that 

J PidiP) = j PidiP) j P(dx)e^^^'^^ (3.31) 

and if, if < X >o= / P^{dip)X, it holds 

< V'x,iV'+i >o= I < ii^i'^ + iV'4'^)(4'^ - i^^^) >o=< >o (3.32) 

as < ip^^^ip^'^'^ >o= if A = (the two systems are decoupled). In the same way 

< V^x.-iV'^-i >o= I < (V^4'^ + *^)(V5i'^ - ii^P) >o=< i^L'^i'^ >o (3.33) 
and finally 

< V'x,iV'J,-i >o= I < (V'i'^ + i^^^W^'^ - #f ) >o=< ^i'^i'^ >o (3.34) 

Note that < ip^ipy >o are the Ising correlation computed in [MPW], [SP] where 
is found that 

< >o= ^ E ^^^^ci,i(k) (3.35) 

k 

< >o= ^2 E ^^^^^2,2 (k) (3.36) 

k 

< i^tW >o= ^ Yl ^^^^ci,,(k) (3.37) 



M2 

k 



where 



A = 2t(l - t'^){2 - cos ko - cos /c) + (t + 1 + V2f{t - ^2 + 1)^ (3.38) 

and 

ci,i (k) = ^ [(2t + 2\/2t) (sin ko cos A; — i sin A; cos ko) + (2\/2 + 3 + 1^) (sin ko — i sin k)] 

(3.39) 
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C2,2(k) = ^ [{2t+2V2t) (sin ko cos k-i sin k cos ko) + {2V2+3+t'^){- sin fco -i sin k)] 

(3.40) 

ci,2(k) = -ci,2(k) = ^[(2V2t^ +4t^) cos /coos /co+ 
if + 2tv^ + t) (cos /c + cos /co) - 2 - 2v^ + 2t2] (3.41) 

Note that 

CM(k)c2,2(k) - C,,2(k)c2,l(k) _ ^ , ^ ^)2 ^ ^^^^^ 

A(k) 

where /(k) = O(k^) and with bounded derivatives. 

4. Renormalization Group for light fermions 

4.1 Multiscale analysis. We start considering Zl^j~~~ in the sum (2.22), written 
as in (3.8). We introduce a scaling parameter 7 > 1 and a positive function 
X(k) e such that 

x(k)^x(-.)^{j (4.1) 

where 

|k| = y/sin A;^ + sinA;2 . (4.2) 
We define also, for any integer h <0, 

A(k)=x(7-'k)-x(7-'^+'k); (4.3) 

we have, for any Hm < 0, 



X(k)= A(k) + x(7-'^"k). (4.4) 

h=hM+'i- 

Note that, if /i < 0, /^(k) = for |k| < ^07^""^ or |k| > to7''+S and /^(k) = 1, if 
|k| = to'y'^. Therefore 

fhO^) = yh < Hm = mm{h : toi^+^ > V2{TrM-^)^} , (4.5) 

and 

1 

1= J2 fhi^) /i = l-x(k). (4.6) 

h=hM 
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We define a sequence of effective potentials V^^^ip) defined iteratively in the 
following way; assuming that we have integrated the scales h= 1,0,-1,-2,..., h+1 



^-M'^Em 



J i'^.,.n„c.(#(^'^))e-^"^^^^^^'''^)-^^^'' , V('^)(0) = , (4.7) 



where 

Pz„m„cAdi^^- 0= 11 11 JfT^. 

k:C-Vk)>0'^=±l 

(4. 

t k:C-^(k)>0 u>,uj'=±l 

aik)-' = J] fjik') , (4.9) 

j=hM 

and the 2x2 matrix T/i(k') is given by 

1 fZ{ismk + smko)+fi,i{k)Z-^ -imft_i(k) + /i,2(k)Z-i 

Co + /o,o(k) V imh-i{k) + f2,i{k)Z-^ Z {i sin k - sin ko) + f 2,2(^)2 

(4.10) 

with nil = Co{t — tc) and Co, Z defined in Theorem 2. 
Finally V^'^) is given by 

cxD * 2n 

vW(v,(^'^))= J] n^"^^S"^2iL(^i'-'^^n) (4-11) 

4.2 r/ie localization operator. We define an £ operation, for h < 0, in the 
following way: 

1) If 2n = 4, then 

>^^<a,c,(ki, k2, ka) = ^++' ^++) ' (4.12) 

where 



2) If 2n = 2 then 



4 

r7,r7'=±l 



l//ebbroio/2008; 22:32 



27 



• +M'!^,Jk,,^)iv-sink + rj'^smko)] (4.14) 
3) In all the other cases 

CW^n,aJ^U---Mn-l) = 0. (4.15) 

By (4.12) the operator C satisfies the relation TIC = 0. 

Remark. First note that in the limit M ^ oo (4.14) becomes simply 

^W^'i,M = r£U(0) + kodkM'ijO) + kdkW^^^lm, (4.16) 

hence £M^2'^^^(k) has to be understood as a discrete version of the Taylor ex- 
pansion up to order 1. Moreover the localization operator acting on the effective 
potential in the x-space representation can be seen as an operator on monomials 
of Grassmann variables; for instance (4.12) implies in the x-space representation, 
in the M — > oo limit 

^ E tli^^^7Srwt\^i,...,^,) = W40,0,^ (4.17) 

xi,...,X4i=l i=l 

where the r.h.s. of the above equation is always vanishing unless Y[i=i V'4f,j'f * is a 
permutation of 'i/'if,+^^if,+ '^^^-'^ '■> hence JC acts on a quartic monomial 
producing a local expression. Analogous considerations can be done for n = 1. 

4.3 Wc have, before continuing, to exploit the consequences of some symmetries. 
There are no local terms of the form ^/'x,i'*/'x,i) such local terms can be written as 
V'x^Vx^'*) but the model is invariant under the transformation 

^(1), ^(1), ^(1), ^(1) ^ -x^^), -x(^) 

^(2), ^(2)^ ^(2)^ -(2) _ ^^(2)^ ^-(2)^ ^(2)^ -(2)^ ^4_lg) 

hence such terms cannot be present as they violate such symmetry. 
There are no local terms of the form ';/'i,x'*/'-i,x or '*Ai^xV'ii,x (^^ '0i,xf^'0-i,x) 
V'tx^V'li.x) ; in fact, 

V'i,xV'-i,x = li^l^^^^^^^ - ^I^^^W + i^^^W + i^^^W] (4.19) 

and the last two terms violates the symmetry (4.18); moreover the first two terms 
are odd in the exchange (1), (2) — > (2), (1) and the model is invariant in the ex- 
change (1), (2) — i> (2), (1). Moreover the model is invariant under complex conju- 
gation and the exchange 

i;^^ , i;^^ ^ i;^^ , i;^^ x^"^ , X^"^ ^ xi"^ , X^"^ ; (4.20) 
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this follows from the fact that, from (2.31), H^°'\H^'^\ V^°'\ V^°'\ written in terms 
,X^^\x^^\ arc invariant under such transformation. Hence the co- 
efficient of the local part of the quartic (non vanishing) terms is real; in fact 

t()(0, 0, O)'0J*^xV'i,xV'-i,x''/'-i,x = 'fi'(0, 0, 0)'i/'4^'*'i/'4^Vx^'*V'x^'' must be equal, by the 
above invariance, to u'*(0, 0, 0)V':L^Vx ^V'x Vx \ hence u)(0,0,0) = zD*(0,0,0). Fi- 
nally the combination of local terms "iAx iV'x -i+'*/'x -iV'x i equal to ^[ip^'^^^^ — 
V'x^''^x^'*] so it cannot be present as it violates the symmetry (4.18). On the other 
hand 'i/'x.iV'x,-! ~ V'x,-iV'x,i is equal to 'il)^'^ 'i/'x^''V'4^'*]; hence the coeffi- 

cient of the local part is imaginary; in fact w{Q)[%Ij'^^ ■ip'^'' + V'x^''V'x^''] must be equal 
to w*(0) ['i/'i^Vx^'' + '0x^'*'0x^''], by the invariance under complex conjugation and 
(4.20), hence w(0) = -w*(0). 

4.4 By (4. 12), (4. 14), (4. 15) and the symmetry relations in §4.3, we can write 
CV^^^ in the following way: 

where, if |A|, |z/| < £, -Sj^t^^i, z^^^^i = 0(£), Si = 0(£), h = 2AsechV(a + 6) + 0(£2), 
= u + 0{e) and 

F{<h) _ V — V ib^-^^+ib^-^^- 

~ ^k,u^ S^k-u^ ' 

a;=±l keX>ju 

^ ^ ki,...,k4ex'M 

kex»M kex»M 

^;-''=F E /.we^'^'i". ^r''= E f^*^'-^ 

kex»M kex>M 

(4.22) 

with /i = sin/c and /2 = sinfco- The constants Hh, Shjlh, Zh,C(h are real. At the 
end of our iterative construction it will appear that it is possible to write (see 
Remark 1 in §4.9) 

W^ilU = ^2ll + <£l (4.23) 

with W2^^^^ vanishing if at least one ruk = 0, for 1 > k > h + 1, and W2^^^^ is the 
rest; we define 



(4.24) 



?7,77'=±1 77,77'=±1 
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In the same way we include in Yl'i=i Z^u- ^i,u>,hFi~^^ all the terms not in ZhF^-^^ + 
a^F^^^ with Zh,,afi real. 

We renormalize the free integration Pzh,mh,Ch{^'^^~^^) by adding to it part of 
the r.h.s. of (4.21). We get 

(4.25) 

where mh-i Ch^^'^^~'^^'^ obtained from Pzh,mh,Ch{d'^^-^^) by substituting 

Zh with 

Zh-i{k) = Zh[l + Cj;\k)Z-\Co + fo,om^h] (4.26) 

and m;i(k) with 

m„-i(k) = ^^K(k) + Cj^\k)iCo + /o,o(k)).;.] ; (4.27) 

moreover V^^^ = V'^^^ - ZhShF^-^'^ - ZhZh{F[-^^ + F^^'^). We will call mfe(O) = 
mh- The r.h.s of (4.25) can be written as 



(4.28) 

where the factor exp(— M^t^) in (4.25) takes into account the different normaliza- 
tion of the two integrations and 

= Z^(l + CoZ-hf,) , A(k) = Z^_i[^i^-4 - ^i^] . (4.29) 
Note that //i(k) has the same support of //i(k). The single scale propagator is 

/ i-..-.,....,/,' (#"*') = . (4.30) 

where 

- y) = Ee"^'^^^"^^A(k)[r."Hk)].,.' , (4.31) 

k 

and T^^(k) is the inverse of the T^(k) defined in (4.10). The large distance 
behaviour of 5'([f^/(x — y) is given by, if \zh\ < \ and sup^>/j \ -^^\ < e''^', given 
the positive integers uq, ni and putting n — hq + rii, and the constant Cjv^„ 

- y)i < ^^^- i + (^/.^d(x-y))r ^'-''^ 
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where dx denotes the discrete derivative. It will be useful to write 

gL'li^ - y) = 9i%J^ - y) + - y) + aL'U^ - y) (4-34) 

with 

giU,.(x-y) = ^E^-'"'-'" . .""fi". . ("5) 

-'^^ ^ —Zusmko + tZsmk 

obeying to the bound (4.32) while 

\SZBT&^ - y)l < gA'." i + (-,4d(x-y))r 

- y)| < CM.nl^l'^^pp^g^^^^ . (4.37) 



Moreover 

9i%i^ - y) = ^'^(x - y) + Jx - y) (4.38) 

with 



-'w! Z ^ sm /Co + -Z^ sm /c^ + (k) 

verifying (4.33) and ^^'^^^^(x — y) verifying (4.36). 
We now resettle the field so that 

v('^)(V^V^-'^^) = v^'^Ha/^V'^-'^^) ; (4.40) 

it follows that 



.(<h) 

(4.41) 

where 



i^h = -^^rih , Sh = -^^{ah - Zh) , Xh = i^^^flh , (4.42) 

and Ci,a;,/i = ■^Zi,u>,h, Ci,u>,h = We CaU Vh = {Xh, Sh, \Ji,^ Ci,u>,h, U,a. Ci,u.,h) 

and the set of {vh,mh,nh) are called running coupling constants. 
If we now define 

^-v(^-^HV^^'-''-'')-^f3^r. = [p ,_l(dW,W)e-^^'H^/^^^-''^) , 

(4.43) 
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it is easy to see that V^^-^\^/Zh-i^'^-^~^^) is of the form (4.11) and that 

Eh-i=Eh + th + Eh. (4.44) 
It is suflBcient to use the well known identity 

oo ^ 

M^^, + V('^-^)(V^V'(^'-^)) = E A(-l)"^'^r''(^^'^(A/^V'^^'^))) , 

n=l 

(4.45) 

where EJ^"^ denotes the truncated expectation of order n with propagator Z^^^g^^^^,, 
see (4.30), and observe that ^^-''^ = ^(<'^-i) + ^C^). 

Let us define 

h* = inf{/i : > /t > /iM, ^07''"^ > 4|m;,|,V/i :0>h>h}. (4.46) 

Of course this definition is meaningful only for mo such that \mo\ < 
The integration of the scales from h* to Hm will be performed "in a single step" . 
This follows from thw following result 

4.5 Lemma. Assume that h* is finite uniformly in M , so that |m/i*_i7~'* | > R, 
for a suitable constant R and define 



Then, given the positive integers N, no, ui and putting n = uq + ui, there exist a 
constant CN,n such that 

Kd^'9l-::P{^; y)l < giv,n ^^(^,?j^(7-y)|)^ • ^^-^'^ 



4.6 Remark. Let us now explain the main motivations of the integration proce- 
dure discussed above. In a Renormalization Group framework one has to identify 
the relevant, marginal and irrelevant effective interactions. By a power count- 
ing argument one sees that the terms bilinear in the fields are relevant, and the 
terms quartic in the fields or quadratic with a derivated (in the x space) field 
are marginal. As it was shown in §4.3, there is only one kind of relevant term, 
and this is a rather crucial point. The Renormalization Group flow of the rele- 
vant terms can be controlled by the introduction of a counterterm; as we have at 
our disposal only one counterterm, it is important to have only one kind of rele- 
vant effective interactions. The unique relevant effective interaction has the form 
V'xjiV'x,-! ~ V'^.-iV'x,!) which can be interpreted as a mass term in a fermionic 
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quantum field theory fike the Thirring model. There is however an important 
difference: in fermionic models coming from QFT if there is no mass term in the 
formal free action then no mass terms are generated by Renormalization Group 
iterations, in absence of spontaneous symmetry breaking; the reason is that such 
models are invariant under the local gauge transformation •i/'x.oj ~^ e*'^""V'x,a; 
there is no mass term in the formal free action. In our model this is not true, 
as the interaction is not invariant under this symmetry; hence even if t = tc (or 
mi = 0) a mass term can be generated in the RG iterations. Hence we collect all 
the relevant terms which are vanishing if rrik = 0, 1 > k > h+1, in s^, which we in- 
clude in the fermionic free integration; the "mass" has a non trivial flow producing 
at the end the critical index of the correlation length. The remaining terms arc left 
in the effective interaction; they are constituting the running coupling constant Uf^ 
(in a gauge invariant theory i/h = 0) whose ffow is controlled by a suitable choice 
of the counterterm v; the running coupling constant i/h takes into account the 
renormalization of the critical temperature. The running coupling constant Vh are 
dimensionally marginal; again the Ci,w,/i7 Ci,oj,h correspond to effective interactions 
which are absent from the free action, but which could be possibly generated by 
the RG iterations, and they would be absent in a QFT fermionic model. Finally, 
due to the mass gap, the propagator of the integration of all the scales between h* 
and Hm has the same bound as the propagator of the integration of a single scale 
greater than h*; this property is used to perform the integration of all the scales 
< h* in a single step. 

4.7 Theorem 3. Let h > h* >0 and, for some constants ci, if 



max[|t)fe| + \uk\ 

k>h 



sup 

h'>h 



rrih' 



rrih'-i 



< e 



sup 

h'>h 



Zh' 



< e 



(4.49) 



there exists a constant e (depending on ci) such that, if £h < £, then, for a suitable 
constant cq, independent of ci, as well as of M, the kernels in verify 



|W^f)(xi,...,X2„)| <MV'''''=^''''°nco£ft)"'""^''"-'^ 



(4.50) 



Xi ,...,X2n 



where 



Dk{Pvo) = -2 + n + k 



(4.51) 



and k = YlfZi 



4.8 Proof. We write V*-'^-' in terms of a tree expansion, similar to that described, 
for example, in [BM]. 
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h h + 1 K 0+1+2 

We need some definitions and notations. 

1) Let us consider the family of all trees which can be constructed by joining 
a point r, the root, with an ordered set of n > 1 points, the endpoints of the 
unlabeled tree, so that r is not a branching point, n will be called the order of the 
unlabeled tree and the branching points will be called the non trivial vertices. The 
unlabeled trees are partially ordered from the root to the endpoints in the natural 
way; we shall use the symbol < to denote the partial order. Two unlabeled trees 
are identified if they can be superposed by a suitable continuous deformation, so 
that the endpoints with the same index coincide. It is then easy to see that the 
number of unlabeled trees with n end-points is bounded by 4". We shall consider 
also the labeled trees (to be called simply trees in the following); they are defined 
by associating some labels with the unlabeled trees, as explained in the following 
items. 

2) We associate a label h < with the root and we denote T^^n the corresponding 
set of labeled trees with n endpoints. Moreover, we introduce a family of vertical 
lines, labeled by an an integer taking values in [h,2], and we represent any tree 

n SO that, if V is an endpoint or a non trivial vertex, it is contained in a 
vertical line with index hy > h, to be called the scale of v, while the root is on the 
line with index h. There is the constraint that, if v is an endpoint, hy > h + 1; if 
there is only one end-point its scale must be equal to h + 2, for h < 0. 
The tree will intersect in general the vertical lines in set of points different from 
the root, the endpoints and the non trivial vertices; these points will be called 
trivial vertices. The set of the vertices of r will be the union of the endpoints, the 
trivial vertices and the non trivial vertices. Note that, if vi and are two vertices 
and vi < V2, then hy^ < hy^. 

Moreover, there is only one vertex immediately following the root, which will be 
denoted vq and can not be an endpoint; its scale is /i + 1. 

3) With each endpoint v of scale hy = +2 we associate one of the contributions 
to V*^^-* given by (3.9); with each endpoint v of scale hy < 1 one of the terms in 
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^yC^f-i) defined in (4.41). Moreover, we impose the constraint that, if v is an 
endpoint and hy < 1, hy = hyi + 1, if is the non trivial vertex immediately 
preceding v. 

4) If V is not an endpoint, the cluster Ly with frequency hy is the set of endpoints 
following the vertex v; if v is an endpoint, it is itself a {trivial) cluster. The tree 
provides an organization of endpoints into a hierarchy of clusters. 

5) We introduce a field label f to distinguish the field variables appearing in the 
terms associated with the endpoints as in item 3); the set of field labels associated 
with the endpoint v will be called ly. Analogously, if v is not an endpoint, we shall 
call /„ the set of field labels associated with the endpoints following the vertex v; 
x(/), e{ f) and uj{ f) will denote the space-time point, the e index and the uj index, 
respectively, of the field variable with label /. 

If h < 0, the effective potential can be written in the following way: 

CO 

V^''Hy%:i;^^'''>)+M^Eh+i^Y. E V^^'Hr, ^ni^^^'^^) , (4.52) 
where, if vq is the first vertex of r and ti, .., (s = Sy^) are the subtrees of r with 

root Vq, 

V^'^\t, y/Zhip^-^^) is defined inductively by the relation 



vW(t,v^V^^'^)) = 



and V('^+i)(r„yz;:V'^-^+^^) 

a) is equal to TZV^^+^^Ti, ^/Zi^'ip^^^+^^) if the subtree n is not trivial; 

b) if Ti is trivial and h < —1, it is equal to one of the terms in jCV'^^+^^ {y%^'ip(^^+^^) 
or, if /i = 0, to one of the terms contributing to V^^^ifj-^) (3.9). 

In (4.53) ^"^^1 denotes the truncated expectation with respect to the measure 
UcoP(di)l^^^^). Ther.h.s. of (4.53) can be written more explicitly in the following 
way. Given r G there are n endpoints of scale 2 and only another one vertex, 
Vq, of scale 1; let us call vi, . . . the endpoints. We choose, in any set Iy^, a 
subset Qy^ and we define Py^ — UiQy^. We have 

VW(r,V'(^°)) = 5;vW(r,P„J, (4.54) 

V^'\r,Py,) = J d^yj^'{Py,)Kl]lj^y,) , (4.55) 
1 " 

4X(x.o) = :^^^ni^^'HPv^QvJ,...,i^^'\Pv^QvJ]ll^^fM , (4.56) 

i=l 
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where we use the definitions ^^^^\Py) = n/eP, ^""^-^Vif//,!^//), and i^^?(x„.) 
are the kernels in (3.9). We now write V(°) as £V(°) + 7^V(°^ with £V(°) defined 
as in (4.21), and we write for T^V*^*^-* a decomposition similar to the previous one, 
with nV^^\T,Py^) in place of V(°) (^. ^^'o)- By iterating the previous procedure, 
one gets for V^'^\t,iIj^-^^), for any r G T^^n, the representation described below 
(see [GM] and [BM] for details). 

We associate with any vertex v of the tree a subset Py of 1^ , the external fields 
of V. These subsets must satisfy various constraints. First of all, if v is not an 
endpoint and vi, . . . , Vg^ are the vertices immediately following it, then Py C IJiPvi ', 
if V is an endpoint, Py = ly We shall denote Qy^ the intersection of Py and Py^, 
this definition implies that Py = UiQy-. The subsets Py\Qy-, whose union X„ 
will be made, by definition, of the internal fields of v, have to be non empty, if 
Sy > 1. Given r G there are many possible choices of the subsets Py, v & t, 
compatible with all the constraints; we shall denote Vr the family of all these 
choices and P the elements of Vr- We can write 

VW(r,V,(^'^)) = ^ V('^)(t,P) . (4.57) 
PeVr 

V^'^\t, P) can be represented as 

V^^\t,P) = V^"'"'' y" d^yJ^^'^\Py,)Ki'^+'\^y,) , (4.58) 

with K^'^p^^Xy^) defined inductively (recall that hy^ = /i + 1) by the equation, 
valid for any v & t which is not an endpoint, 

4'P^(X.) = -,{^)'^fl[Ki':^^'\^V.)] 8l[^^^^\Py,\Qy,),...,^^^^\Py^^\QyJ] , 

(4.59) 

(2) 

Moreover, if v is an endpoint, Ky (xy) is defined as in (3.9) if hy = 2, otherwise 

{Aft„-i5(xi - X2)5(x2 - X3)5(x3 - X4) if V is of type A, 

j'^^-'.n.-iSi^-y) if ^ is of type., 

d/i^_id(x - y) ^ if ^; IS of type 5, 

C/i.-i,u.,i<^(x-y), C/i.-i,u.,i^(x-y) if is of type z,^ 

If Vi is not an endpoint, 

Ki':^+'\^v,) = 7^<pt?(x..) , (4.61) 

where is the subtree of r starting from v and passing through Vi (hence with 
root the vertex immediately preceding f ), P*^*^ are the restrictions to Tj of P. 
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(4.57) is not the final form of our expansion, since we further decompose V^'*^ (r, P), 
by using the representation of the truncated expectation (3.21). If we apply the 
expansion (3.21) in each non trivial vertex of r, we get an expression of the form 

V('^)(r,P) = ^^''^'^ J2 [ ^x„,V;^-'n^'.o)<'i^,T(x.o) , (4.62) 

where T is a special family of graphs on the set of points x^;^ , obtained by putting 
together an anchored tree graph T„ for each non trivial vertex v. Note that any 
graph T e T becomes a tree graph on x^^,, if one identifies all the points in the 
sets Xy, for any vertex v which is also an endpoint. 

Taking into account the efi'ect of the TZ operation we obtain (see [BM], §3 for a 
detailed proof of a similar formula) 



W^r,P,T(x„J = J2 



aeAx 



n 

hi 

V* 

.i=l 



n 77 f dPrMdetCt-^'-iK) 

V J 



not e.p. 



not e.p. 

(4.63) 

where "e.p." is an abbreviation of "endpoint" and, together with the definitions 
used before, we are using the following ones: 

1) At is a set of indices which allows to distinguish the different terms produced by 
the non trivial TZ operations and the iterative decomposition of the zeros; 

2) t;^, . . . , f* are the endpoints of r and hi = hy* ; 

3) ba{l)-, Qq;(//~) and qa{fi~) are positive integers < 3; — /; 

4) if V is a non trivial vertex (so that Sy > l),the elements GJl-l,^, of G^'''^^(t„) 
are of the form 

G^M'r = tM'^''"^^^^^^'"^^^^^'^^"^^'"^^"^^'"^^i'-,l+('^i.- - Yi'r) ; (4.64) 

if V is trivial, T„ is empty and J cZPt^ (t^) det G^'"-^'' (t„) has to be interpreted as 
1, if |Xt,| =0 {Ty is the set of internal fields of v), otherwise it is the determinant 
of a matrix of the form (4.64) with ti^i' = 1. 

It would be very difficult to give a precise description of the various contributions 
to the sum over At, but fortunately we only need to know some very general 
properties, in particular that \At\ < C" for some constant C and that for any 
a e At, the following inequality is satisfied 

f€lvQ l€T vnot e.p. 
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where ha{f) = hy^ — 1 if / e P^^, otherwise it is the scale of the vertex where the 
field with label / is contracted; ha{l) = hy, if I E Ty and 

ri if|P,| = 4, 

ziPv)=< 2 if |P„| = 2, (4.66) 
1 otherwise. 

By a standard computation (see for instance [BM] , §3) by bounding the determi- 
nant by the Gram-Hadamard inequality (see (3.23)) we obtain 

J d^y,\Wr,T>,TM\ < C^M'elr'''''^''^^^ ■ 

(4.67) 



n {j^c^^- 



V not e.p. 



with — 2 + + z(Py) > 0. In order to perform the sums note that the number 
of unlabeled trees is < 4"^; fixed an unlabeled tree, the number of terms in the 
sum over the various labels of the tree is bounded by C"^, except the sums over 
the scale labels. In order to bound the sums over the scale labels and P we first 
use the inequality 

\Pv\ 



(4.68) 

V vnot e.p. 

where v are the non trivial vertices, and v' is the non trivial vertex immediately 
preceding v or the root. The factors ry-'2a{hy-h^,) r.h.s. of (4.68) allow to 

bound the sums over the scale labels by C"'; a is a suitable constant (one finds 

Finally the sum over P can be bounded by using the following combinatorial 
inequality, trivial for 7 large enough. Let {pv,v G r} a set of integers such that 
Py < '^iliPvi for all V e T which are not endpoints; then 

^7-^<C". (4.69) 
II not e.p. Pv 

It follows that 

E n ^"^^ n E7-"<^^- (4.70) 
P x; not e.p. I! not e.p. Pv 



4.9 Remark 1. The decomposition in (4.23) respects the determinant struc- 
ture of the truncated expectations, as we can decompose the propagators as in 
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(4. 34), (4. 38), obtaining, for any tree with n end-points, a family of C"' different 
contributions to M^2'*^ which can be bounded as before. Hence we can include in 
W^"^^ the terms with at least a propagator g^^_^, defined in (4.39). 
Remark 2. If the tree r has an end-point on scale k, the bound (4.50) can be 
improved by a factor it follows immediately from (4.68). In particular 

if to an end-point is associated a term TZV^^'^ there is an extra factor 7"^^; this is 
a short memory property. 

4.10 The above results were proved for ^2}' ' '"; a similar analysis can be 

repeated for Z|j ''^ , for any value of e^^\e ^^\e^'^\e The only dif- 

ference is that one has in addition the function Qg in the interaction and the 
oscillating functions e*P^^. One can split in a part identical to the one for 
called V'^-'^^ and the rest; hence we repeat the multiscale analysis by 
writing V^'^^ = V^'^^ -\-V with V given by a sum of trees with at least an 
end-point associated to V^-*^) — V^^^; we define the localization operators acting non 
trivially only on V'^'*^ (and defined as above). It is easy to see that the terms from 

trees with at least one end-point associated to V*-^-* — V*^^-* are vanishing in the 

-h* 

limit M — 00; in fact the bounds for such terms is improved by the factor ^ , 
by simply dimensional considerations, and we will see in the following section that 
7~'* is a finite number independent from M, hence such terms are vanishing in 
the limit M — > 00. 

5. The flow of the running coupling constants 



5.1 By the analysis of the preceding section it follows that the running coupling 
constants Vk.i'k^'n^k^ 1 > k > h* , verify a set of recursive equations called Beta 
function equations, of the form 

Uh-i = 'ji^h + (^y{vh, T^h\ vi, Vl) 

Vh-i = Vh + (3y{vh, i^h; ^^1, t^i) (5.1) 



rrih-i 
rrih 



Zh-l 1 , ohr \ 

By repeating the analysis for proving Theorem 3 to the functions we have that, 

if (4.49) holds, then the P'j^ are expressed by convergent series. We want to show 

_ 1 

that there exists a positive constant as such that, if |t — td > e , it is possible 
to find a function ui (hence a function u) so that (4.49) holds (if A is small enough). 
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Iterating the first of (5.1) we find 



k=h+l 



where now the functions are thought as functions of z^^, . . . , z^i only. 
If we put h — h* in (5.2), we get the foUowing identity: 



i^h* = 1 



-h* + l 



k=h* + l 



and we look for a vi verifying 



(5.2) 



(5.3) 



so that 



k=h* + l 



h 

^h = -l-^ E 7^-V^(z.fc,...,i.i), h*<h<l 

k=h*+l 



(5.4) 



(5.5) 



Let he u = {i>k,h* < h < 1} and ||z/|| = supf^*^f^^i \i>h\; we call M the set of 
i> with bounded norm. We decouple the beta function equations (5.1) imagining 
that in the last three equations of (5.1) Uk is replaced by acting as a parameter; 
we call mh{i>),Vk{i^) the solution of the second and third of (5.1) as functions of 
the parameter u. We shall prove the following lemma. 



5.2 Lemma. There exists e and such that, for \t — tc\ > e°-^^i , |A| < £ and any 
V such that < Cj^|A|, it holds, if mo — \t — tc\ and C,C^, are constants 



|Aft(z>) - All < ^ \Xh{i>) - Xh{i>')\ < C\X\ max - 

4 k>h 



(5.6) 



\Sh{i>)\<^-^ Mi>)-5h{i>')\<C\X\ 



max li/k — i^kl 

k>h 



\Ci,uj,h{i^)\ < |Co,a;,ft| + |A| \Ci,uj,h{i^) - Ci,u:,h{i^')\ < C\X\ max\uk - u'^ 



k>h 



\Ci,uj,h{i>)\ < |Co,a;,/i| + |A| \C,i,uj,h{i>) - Ci,uj,h{i'')\ < C | A | max | - 



k>h 



--l|<max|^,-^. 



^h-i(z>') 



Zh{i>) Zh{i>') j>h 



< max — 



k\ 
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with a2, a4, ci, C2, C3, C4 positive constants. 

Remark. An obvious corollary of the above statement is that there exits a finite 
(uniformly in M) h* such that mh*-i'~i^ ~^ > « for a suitable constant k. 

5.3 Proof. The proof is done by induction. By iterating the second equation of 
(5.1) we find 

1 

A/,_i(z>) - Ai = ^f3l{vh{i>),i'h\-\v^{i)),v^) (5.7) 

k=h 

and by induction 

\PI{vh{v), vw. ...-My). yx)\ < C7i[|Ai|2 + A?^ + Xlr% (5.8) 

The first addend of the r.h.s. of (5.8) is a bound on the terms containing only 
propagators gi^l^i^ — y), see (4.35), and trees with end-points v with hy < 1, 
and the bound follows from the fact that the second order terms cancels out. 
The second addend is a bound on the terms containing at least a propagator 
guj[^,L02i^ ~^ y) (4. 39), (4. 37) (the boimd follows from the short memory property 
and (4.33), (4.37)). The third addend is a bound on the sum of terms with at 
least a propagator g!^i\uj2 (x — y) or from trees with at least an end-point v with 
hv = 2, and we have used (4.36). 
Inserting (5.8) in (5.7) we get 

|A,-i(^) - All < C^[\h\\X^\' + A? ^^[^ + < ^ (5-9) 

k=h ^ 

for 03 large enough and A small enough, where we have used that Ylt=h "^^^ 
bounded (by induction and the definition of h*). 
Moreover the expansion for A/i_i(z/) — A/i_i(z/') is given by a sum of terms similar 
to the ones for Xh-i{i^) in which a Vk is replaced by ffc(z^) — Vk{i'') bounded by 
CA^ maxj>fc \i>j — i>j \ or or replaced by the relative difference; one finds, for 
as large and A small enough 

|A„-i(i>)-A^_i(i>')| <C2[|/i||A|3 

I I 

+ |A| y][^| + 7'''']] max |z>fc - i>^| < C|A| max \i>k - i>k\, (5.10) 

^ — ' 7*^ k>h k>h 

k=h 

where the first term is a bound on the terms containing only propagators ^i'*^(x — 
y), and trees with end-points v with hy < 1, containing at least a Vk{u) — Vk{i'') 
or a z/fc — z/^ (note also that the first term with end-point v is A^z^^). A similar 
analysis can be repeated for 5^. 
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The recursive equation for Ci,uj,h is 

Ci,o.,h-iii>) = Ci,o.A^) + Ci,o.A>^las + P,) + (5.11) 

with \/3^\ < Cs\X\^ and < Qb"'^ + ^JA^. In Xl{a2 + (3^) in (5.11) are the 
contributions from terms with only diagonal propagators Qi^lji^) and from trees 
with endpoints v with hy < 1; there is necessarily at least a vertex C. In fact 
assume that a contribution has nx end-points of type A and n end-points of type 
v; there are then tt-a — | + f propagators u> = 1 and — | + f propagators 
a; = — 1, so that n must be odd; the total number of propagators 2nx — 1 -|- n is an 
even number; the derivated integrand is odd under the exchange k — > — k and so 
it is vanishing. Moreover there is no contribution of order XhZh (they corresponds 
to tadpole graphs , whose derivative is zero). In R^-^^ are the trees with at least 

(k) 

an end point v with — 2, or with a propagator goj-/,u;2: or with a propagator 

(k) I ~ 

9W,oj2- It holds by iterating, for 03 large and A small enough, |Ci,w,/i| ^ piju-.o! + |A|. 
Moreover 

I I 

\Ci,u,,h-i{i>) - 4c.,/^-l(^>')l < C3 ma^ \vk - i''k\W\h\ + |A| ^^[^ + 7"']] < 

3>h f— f 7'*' 

k=h 

C|A| max jf/fc — I'fcl (5.12) 

j>h 

Similar computations can be repeated for Ci,oj,h- 

By definition (see (4.23)) mh-i{p) is given by a sum over terms with at least a 
non diagonal propagator ^(^'^^^^(x — y), hence 

"^^-ff =l + Afc(-a2 + M (5.13) 

where 02 > is a constant and \hm\ < C|Ai|; then (assuming Ai > 0; similar 
computations can be repeated for Ai < 0), for A small enough 



Then 



^— 2a2Ai 



and 

3 



mH-i{i>)\ < |mo|7^^-("-^) ^^ < |mo|7^"^('^-^) (5.16) 
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Finally, for large enough 

< C4log/3 max \vk{v) — Vk{v')\ < majcluk — i^k\\M~^ i^-^"^) 

0<k<h j>h 

Finally the last of (5.6) are found in a similar way, noting that ^Zh{i>)^ = 1 + 
A|(o4 + hg), with \hs\ < C'lAij. 

-1 



5.4 Lemma. There exists e, as andC^, such that, for \t — tc\ > e°3^ and |A| < s, 
there exists a z^i(A) so that maxfc>/j* |z^fe| < Ci^|A|. 

Proof. It is sufficient to look for a fixed point for the operator T : M — > M, and 
T is defined in the following way, if u' = T(i>), see (5.5): 

h 

^'h = Yl 7'-'/5fc(^fc(^),%,...,i^o(i>),i>o), h*<h<l. (5.18) 

k=h*+l 

We want to prove that it is possible to choose the constant Cj^ > 1 in Lemma (5.2) 
so that, if |A| is small enough, the set .F = {i^ G M : ||f>|| < Ci^|A|} is invariant 
under T and T is a contraction on it. This is sufficient to prove the lemma, since 
M is a Banach space, as one can easily show. 

By (5.18) and Theorem 3 



h+j-1 

I 

j=h*+l 



Ci,.|A| + Ec"|A| 



n=2 



(5.19) 



where Ci^u is a constant depending on the first order contribution {i.e. the tadpole). 
So, for a proper C,^, < Cz^|A|. 

We then show that T is a contraction on J^. In fact, given i>2 £ by using 
Theorem 3 and Lemma (5.2), we see that, for A small enough 

k=h*+l 

K(i^l): ^k,l, vi{ui), Z>i,i) - Pk{Vk{i'2), Z>fc,2, • • • , Vi{D2), Z>1,2)] (5.20) 

and 

h 

k=h* + l j>h 
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with Ci, C2 constants, so that T is a contraction. 



5.5 If a = 6 = one has the case treated in [SP], corresponding to two indepen- 
dent Ising models with a quartic or a next to nearest neighbor interaction. In such 
a case the local part of the quartic terms is vanishing by Pauli principle, so that 
if |P„| =4 one can apply "freely" a first order renormalization obtaining an addi- 
tional ^ for any v such that \Py\ = 4, in the bounds. At each step one 
can include all the quadratic running coupling constants in the free integrations; 
their beta function is Ci,h-i = Ci,h + 0{j°''^X) and Ci./i-i = Ci,h + 0{j'^'^X), as the 
beta function is sum over all the trees with end-points necessarily at scale hy = 2 
( contrary to our case, as the local part of the quartic terms is vanishing, and 
there are no end-points associated to the quadratic terms as they are included in 
the free integration). Hence all the quadratic couplings are 0(A) and they do not 
change the scaling properties of the propagator on a single scale; the qualitative 
behaviour of Cy for A = or A 7^ are then the same for temperatures up to tc- 



6. Correlation functions 



6.1 Flow of observables. We consider the following functional integral 
where V^^^ given by (3.8) and B = + with 

S^^)(V, 0) = / rfx0(l)(x)[V^JiVx,-l + V'x.lV'J-l] (6.2) 

B('\i;, 4>) = J cZx</,(2)(x)[V'+,<-i + V'x^V'x,-!] • (6.3) 
After integrating the fields ...'«/;('*+^\ 0> h> h*, we find 

(6.4) 

where Pzh,mh.,Chid''^^~^^) ^ire given by (4.8) and (4.11), respectively, while 

S(h+i) ((^)^ which denotes the sum over all the terms dependent on but indepen- 
dent of the ip field, and B^^\'4>^-^\ 0), which denotes the sum over all the terms 
containing at least one cj) field and two t/^ fields, can be represented in the form 

00 „ m 

Sih+l)(^)=J2 / d^,-.-d^mSS+'\^l,...,^m)[ll^^'''\^i)] (6.5) 
m=l'^ 1=1 
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OO CX3 „ 

m=l n=l a.iv 



2n ■ 

(6.6) 

m 2n 



i=l i=l 

Since the field (f) is equivalent, from the point of view of dimensional considera- 
tions, to two tj) fields, the only terms in the r.h.s. of (6.6) which are not irrelevant 
are those with m = 1 and n = 1, which are marginal. The localization C is defined 
equal to zero except in the following cases 

C J dxdy(izW^^2(x;y,z)0(^)(x)V;^,iVzVi = J dxdydzW^^2(x; y, z)0(^)(x)V^^,iVx,-i 

(6.7) 

L J dxdydziy^2(x;y,z)0(2)(^)^a^^-.^^ J rfxdydz;y^2(x;y,z)0(2)(x)V;-iVx,-i 



Hence 



7(1) 7(2) 



£SW(^(<M, 0) ^ t^F^^^) + ^Fi^'^) , (6.8) 
where Z^^^ and z'^^^ are real numbers, such that Z^^ = Z^^ = 1 and 

F^' = / <ix,A«(x)[^if' vi?!-;- + , (6.9) 

F.K'" = / rfx,A(^)(x)[V'g'"+V'K'J+ + ^if'-^K';-] . (6.10) 
By using the notation of the preceding section 

= e-^''- j Pz,.,,m,..,Cn.Ad^^^'-'^) ■ (6.11) 

/ Zh-i,mh-i,ff^ \ T ) 1 



where B^'^-^\tl:^^'^-^\(p) and S^'^\(p) are then defined by 



(6.12) 



The definitions (6.8) easily imply that 

^(1) ^(1) ^(2) ^(2) 

= ^[1 + aiA;, + /ii] = ^[1 - aiA,, + /is] (6.13) 
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with < CA^, z = 1, 2; proceeding as in the proof of Lemma 5.2 

7(1) v(2) 
^-\iC4,h ^ h ^ ^-Aicsh. ^Aici/i ^ /t ^ ^AiC2/i _ (6.14) 

The fields of scale between h* and /il,^ are integrated in a single step, and it 
follows that 

1 

5(0) = -LI3El,p + 5^(0) = -LI3El,p + 5] ^('^)(0) ; (6.15) 

h=h* 

As a fields is dimensionally analogous to two external fields, we get, proceeding 
as in the proof of Theorem 3 (see also [BM] §5 for a detailed proof of similar 
bounds in a related case) 

and if = x, y and (t>a — (pi, <p2 

'a0("i)(xi)a0(«2)(x2)a0(«3)(x3)"'^-°-^ ^ ^Al Jl+(^/.|rf(x-y)|)iV 

(6.18) 



(xi)a0(«2) (X2)a0(«3) (X3)a0«4 (^4) 



=0< 



6.2 Correlation functions. The correlation function < a'^'' a'^^ >t, 
where x' = (a; + 1, xq) or (x, + 1) is given by 

If 0(")(x) = cri^2o(T^'l\ ^ each derivative produces a factor 



sech' JH^Hi%^^^ + tanh J + — V (6.21) 

OJv.X + l.Xr, 



'x.;x+l,Xo 
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with V given by (2.25). We define, using (2.30), (2.25) 

<0(«)(x)0(«)(y) >A;,(i),,'(i),,(.),,'(.) = 



Moreover we call 



(6.22) 



< C»(«)(x)C»(")(y) >A,T;sa),s'a),s(^),s'(^)=< 0(«)(x)C»(«)(y) >a;,(i),,'(i),,(2),,'(.) 

- < 0(")(X) >A.,(i),,'(i),,(2),,'(.)< 0("Hx) >A;sa),s'a),e(^),e'(^) (6-23) 

so that 

< 0(")(X)0(")(X) >A,T= J] (-l)'^e(l),e'(l) ^ (_1)'5^(2).,'(2) 

e(i),£'(i) e(2),e'(2) 

z;"'''"' < 0'°'(x)0'°'(y) >A,T,,-...-.., 

Suppose that x and y are fixed to an M independent value; then 

lim < 0(«)(x)0(«)(y) >A,T;e^,en,s-,s'- - < O^") (x)0(«) (y) >A,T;-,-,-,- = 0. 
M — >oo 

(6.25) 

In fact, as explained in §4.10, the l.h.s. of (6.25) can be written as a sum of trees, 
and in each of them there is a factor e^P'^ — 1 or an end-point associated to 
(5ei,g'i,g2 g/2. With respect to the previous bounds, simply dimensional analysis 

says there is now an extra factor in the bounds so it is vanishing in the limit 
M — > oo, if X, y and t — tc are fixed to an M-independent value. We can then 
simply study < 0(")(x)0(")(y) >a,T;- - - - which is given by the Grassmann 
integral (6.22) with 0^°^^ given by (6.21). By performing the change of variables 
(2.31) and (2.59), we get a sum of averages of monomials in the ip and x fields; 
we integrate the x-fields, as discussed in §3, and we obtain a sum of Grassmann 
integrals of monomials in the ip fields; remembering that 

^a^a)^a^(a) ^ ^[^-^i;Z,,^^ty^+^,y + i^^yi^Z, ,y^tj^\^+ (6.26) 

we find the first two of (1.13) and (1.15) by (6.1), (6.13), (6.15) and (6.16), (6.17). 

we include to contributions to the correlation function corresponding to 



l//ebbroio/2008; 22:32 



47 



monomials with six or more fields, and from (6.18), (6.19) the last of (1.13) 
follows. 

Finally the specific heat is obtained in a similar way, noting that we have to sum 
over X— y, and this produces an extra 7^^ in the r.h.s. of (6. 16), (6. 17), (6. 18), (6. 19). 
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Fig.4: Closed polygon in the high temperature expansion 
of the Ising model partition function. 




Fig.3. The Eight-vertex model in the rotated frame, 
with spins relabelled; the solid Unes represent 
the quartic interaction. 
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Fig.2. The spin representation of the Eight-vertex model; 
the soUd points have xO+x even and the void points, 
odd; the sohd lines represent the quartic interaction. 
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Fig.l: The spin representation of the Ashkin Teller model ; 
the dotted hnes represent the Ising interactions, the solid 
lines the quartic interaction. 
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Fig.5. The second lattice and its polygon configurations. 
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Fig. 6. Local identification of Ising configurations and configurations 
over the second lattice. 



